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1. Introduction 

The investigation of atomic bosons with short-range repulsive interactions in a periodic 
potential by using the Bose-Hubbard model has revealed a quantum phase transition 
between two distinct phases: a superfluid and a Mott-insulator, that exists at sufficiently 
low temperatures [1]. The formalism of the Bose-Hubbard model was successfully 
mapped onto a system of cold bosonic atoms in an optical lattice [2|. The superfluid to 
Mott insulator phase transition was experimentally realized |3] and further examined 
and theoretically digested jl]. Continued progresses have focused on systems of 
multi-component Bose- Einstein condensates (BECs) in an optical lattice |5j, where 
diverse topics such as quantum phase transitions of spin-2 bosons [6], two-component 
condensates [7], and spin-1 bosons with coupled ground states [8] are studied. 

An interesting feature characterizing a variety of lattice models mapped onto atomic 
gases is quantum entanglement. Additionally, cold atom based lattice models have been 
identified as ideal candidates for universal quantum emulation of strongly interacting 
many body systems. While a complete understanding of quantum entanglement and 
correlations in an atomic lattice model remains a significant challenge even in theoretical 
terms [Hj, much has been understood for an important type of correlation, the so-called 
spin squeezing, or pseudo spin squeezing. For those systems that undergo quantum 
phase transitions, the presence and the measure of entanglement is important not only 
at the transition point, but also for the different phases of the system. These systems 
show various behaviors, entanglement and disentanglement, coherent and squeezed spin 
states, mode and particle entanglement for different phases that can be controlled by 
interaction types and strengths as well as lattice configurations. 

Squeezed spin states are states whose spin fluctuation in one of the transverse spin 
components is below the standard quantum limit. It was shown in Ref. [19] a spin-s 
squeezed spin state is a correlated state consisting of 2s spin-1/2 particles. This implies 
a potential connection between spin squeezing and entanglement, due to the existence of 
correlations affecting the separability of a system with many spin-1/2 particles [20]. Spin 
squeezing can occur in many models with a variety of atom-atom interactions [211 122] , 
for atomic condensates inside external traps [20], and for atoms inside optical lattices 
[23]. 

In this work, we are interested in the possibility and the condition for spin squeezing 
in the pseudo-spin of coupled ground states in an optical lattice model with spin-1 
bosons. We hope to explore spin squeezing properties of the system carefully studied in 
Ref. [8]. This paper is organized as follows. In Sec. O we review the model system [81124] 
and describe the mapped Bose-Hubbard Hamiltonian in the mean-field approximation. 
The measure of spin squeezing and quantum entanglement that we employ is introduced 
in Sec. 131 The results of spin squeezing for different interaction regimes are presented 
in Sec. HI Finally, we conclude in Sec. [51 
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2. Model System 

The system we study consists of neutral bosonic atoms with hyperfine spin F = 1 in 
an optical lattice. The optical lattice results from the ac Stark shifts of standing wave 
laser fields, which are dipole coupled to atomic electronic transitions. The off-resonant 
coupling induces virtual transitions to electronic excited states, which upon adiabatic 
elimination give rise to level shifts (ac Stark shifts) in the ground state manifold. These 
shifts are proportional to the intensity distribution of the laser light. Additionally 
two-photon Raman like transitions can couple any two Zeeman states within the spin- 
1 ground state manifold, subject to appropriate polarization selections. In a lattice 
of ac Stark shifts from standing waves, the periodic level shift gives rise to band 
structures. When the lasers are linearly polarized, the Zeeman ground state manifold of 
{Mp = —1, 0, -|-1) remains degenerate in the lattice. For more general cases of coupling 
referred to as the A or \^ scheme with suitable polarizations, two alternate ground states 
become coupled and will be denoted as the electronic modes with a = and a = A [8]. 

We assume that atoms will remain in the lowest Bloch bands as a result of 
the relatively large band gap in comparison to their kinetic energies. Within this 
approximation, the atomic field operator can be expanded in terms of the site localized 
Wannier basis. As carefully presented in Ref. [8], we arrive at the model Hamiltonian 
defined on a ID optical lattice as given below, 

Hbh = - X! I^'^l X! '^ii^'^j 

<7=0,A 

(T=0,A i i 

i 

-5j2noi-iJ, J2 Zl^^i' (1) 

i (T=0,A i 

where is the tunnelling parameter, Ua, K, and P are parameters from the repulsive 
density-density interaction of condensed atoms and the spin-exchange interaction. 5 
parameterizes the energy difference between the electronic internal states a = and 
a = A. /i is the chemical potential, dl^^ and dcrj are respectively creation and annihilation 
operators of an atom in mode a at lattice site i and n^i = a^iCiai- 

As discussed in Ref. [8], the various parameters of the above Hamiltonian ([1]) can 
be given in terms of Wannier spinors, and thus they depend on 9, the angle between 
the polarization vectors of the two counter-propagating linearly polarized laser beams 
in the lin-^-lin configuration of an optical lattice. 

In the mean-field approximation ^25j with ip^ = (d^j) assumed real |8J, we substitute 

dlid^j ^ -^aiaaj + aii) - i/^l, (2) 
into the Hamiltonian ([T]), and arrive at 

H^H— 2 E Mia. + at - +E ^^n.- 1) 

(T=0,A o-=0,A 
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+ KfiofiA —{alalaACiA + a\a\aoao) 

-6ho- (3) 

o-=0,A 

a system of many independent sites. In the above, we have omitted the site index i so 
that effectively, the optical lattice model is reduced to a collection of single site problems. 

The basic idea of the mean-field theory (MFT) is to replace the fluctuating exchange 
field by an effective average field in an interacting many-body system. MFT has been 
found not quite reliable to describe critical phenomena especially at low dimensions [9]. 
In the MFT, one ignores the long range fluctuations of the order parameter which 
causes serious errors at the critical points where the fluctuations dominate the mean 
value [To]. Despite these facts, optical lattices have been extensively studied under a 
MFT approach [21 [TT] . The interaction term in the Bose- Hubbard model for the optical 
lattices, e.g. the interaction terms in ([1]), is due to atom-atom collisions which can 
happen only locally, so that it is an on-site interaction. The sole non-local interaction is 
the hopping term, due to tunneling of the atom between the sites. As in the case of our 
spin-1 model, MFT treats the spin-spin interactions exactly while the kinetic coupling 
is treated approximately. 

MFT, as it is used here, based upon the Bogoluibov symmetry breaking background 
field theory. Bogoluibov theory is extended to describe Mott transition by a specific 
decorrelation approximation in a consistent MFT [12]. It can be systematically improved 
by considering bigger clusters (2 sites or more) to employ MFT [13j. Away from phase 
boundaries such an improvement is not essential for us. The fluctuations are due to 
collective excitations of the system. Focusing at zero temperature, and staying away 
from the phase boundaries, one can expect the predicted MFT ground states are well 
estabhshed, since the collective excitations and associated fluctuations would be weaker 
in comparison to the mean-field order parameter. In our investigations we assume the 
reported ground states [8] describe the system in deep quantum phases away from the 
phase boundaries. 

Similar approach, as is done here, to determine the ground states have been 
employed in a more general system that includes external magnetic field as well |14j . 
MFT cannot be used to examine spin-spin correlations among different sites for which 
effective models can be used [I5]. On site spin fluctuations however can be examined in 
MFT to reveal any particle entanglement associated with the reported ground states [8]. 
The question we address here is how the type and amount of the entanglement among the 
particles in a single lattice site would change when the whole lattice system undergoes 
quantum phase transitions and the use of MFT is sufficient for this question. 

Beyond zero-temperature, a generalization of the method is given in Ref. [16]. At 
non-zero temperatures it is more crucial to test predictions of MFT for low dimensional 
systems against numerical tests. For spin-1 systems, detailed numerical studies became 
only very recently available [17]; but they have ensured that similar level of agreement 
between the MFT predictions and numerical studies as in spinless systems do occur for 
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the case of spin-1 systems. 

In order to test the vahdity of MFT that we use in our model, we studied a simple 
lattice model having two sites. We used the Bose-Hubbard Hamiltonian in ([T]) and i 
runs from 1 to 2 with the periodic boundary conditions. The purpose of this calculation 
is to investigate the effect of inter-site interaction on the single-site state. The exact 
ground state calculations were done by using those parameter values corresponding to 
n = 1 and n = 2 Mott phases in the phase diagrams both for the ferromagnetic and 
antiferromagnetic regimes in the case of ^ = and for a small 6 value. Once the exact 
two-site ground state is determined, we calculate the one-site density matrix by tracing 
out the other site. Following this procedure, the overlap of ground states from MFT 
and exact two-site model can be computed. Our results show that most of these overlap 
values are above 0.95, confirming the success of MFT in calculating one-site ground 
states and so that using it to quantify correlations among particles in a single site. 

In general, many-body wave functions are too complicated to express explicitly, 
but MFT allows for writing down analytical wave functions of the ground states and 
hence one can gain valuable insights into the quantum correlations in such complex 
many-body systems such as spinor condensates in optical lattices. This insight should 
serve as a guide even for comprehending quantum correlations among the lattice sites 
which require beyond MFT calculations, but can still be performed through perturbative 
examinations of mean field ground states. We hope to investigate this in near future. 

A general spin-1 system is described by the symmetry group SU(3). In the model 
considered here, a reduced two-mode description for the two coupled ground states is 
represented by a pseudo-spin- 1/2 algebra, effectively the isospin subgroup of SU(3) [8]. 
The corresponding generators of the SU(2) isospin algebra are given by ^ 

Ti = ^{a^aQ + aJ^A), 

f3 = ^{aldo-a{aA), (4) 
in terms of which the mean-field Hamiltonian ([3]) can be expressed as 

a=0,A 

+ + {K- |P|)f 2 + {K+ |P|)f| 

H n — h At H V - \ n 

8 V 2 4 2; 

_(^^,)f3^^„f„ (5) 

where Af/ = Uq — U\, Us = Uq + Ua, and n = uq + n\. Spin dependent interaction 
terms in this Hamiltonian emulates that of the generalized Lipkin-Meshkov-Glick (LMG) 
model [271 126], or its special case of the two-axis twisting model |19j. Such models are 
capable of generating spin squeezing [I9] and multiparticle entanglement [21], [26] . Our 
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model above, includes tunnelling and collision effects in addition to the generalized LMG 
interaction terms. 

When the lattice parameter ^ = 0, the two modes have the same energy and 
Jo = Ja = J , Uq = U\ = U , K = U + P , and 6 = [8j. The simplified Hamiltonian ([3]) 
takes the following form 

7^af(f) = - 2 J E [{ai + a^)^^ - 

cr=0,A 

+ 2 (f/f 2 + Pf|(y) + ah, (6) 

for both antiferromagnetic (P > 0) and ferromagnetic (P < 0) interactions [24j, where 
we have used = + h/2 for the collision interaction in terms of the total isospin 
operator with a = —3U /2 — P/2 — /i. The spin interaction now reduces to that of a 
single-axis twisting type [T9j. 

The above considerations show that our model allows for the investigation of effects 
due to tunnelling and collision on spin squeezing induced by either the two-axis twisting 
interaction as in the generalized LMG model or the single-axis twisting interaction in 
the simplified case. In the general case of the LMG model, particle entanglement thus 
exists for atoms in the non-degenerate ground state modes, which become degenerate 
for the special case of a lattice with ^ = 0. 



3. Spin Squeezing and Quantum Entanglement 

Squeezed spin states defined by Kitagawa and Ueda [19] is widely used in atomic physics, 
especially in the context of particle correlation and entanglement. A criterion was found 
recently connecting many atom entanglement and correlation originally from atoms in 
a Bose-Einstein condensate (BEG) [20]. If the squeezing parameter 

is smaller than 1, the two mode bosonic many atom state under consideration is spin 
squeezed along the direction of a. J is the total pseudo spin operator, while a, /3, 
and 7 denote three orthogonal axes. The condition for < 1 coincides with the non- 
separability criterion of a density matrix for two state boson [20]. Thus can be 
used to measure quantum entanglement in the two state atomic system discussed above. 
In our study outlined below, we examine spin squeezing for the on-site isospin algebra 
by calculating the variance and expectation values of the corresponding generators Tj 
defined in (jlj). Our results show clearly the existence of quantum correlations between 
atoms on the same lattice site. 

To identify pairwise entanglement in our many-body system, we can make use of 
a direct relationship between concurrence [28], which is well-known and represents a 
widely accepted measure of bipartite entanglement, and spin squeezing criterion [29] . 
Thus, we take ([7]) as an indicator for two-particle entanglement. We will in addition 
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also calculate the concurrence and compare the results with the squeezing parameter 

In view of the significant difficulties of measuring spin squeezing along any arbitrary 
direction a, our investigation will focus on the simplest case of a single orthogonal 
configuration with three fixed axes. Other orthogonal axes configurations may be 
sequentially searched for if the optimal squeezing is to be found. For this aim we 
only need to rotate the coordinate system about each of the axes by an angle 0. For 
example if the rotation is about the axis-3, ^| remains the same, while the squeezing 
parameters for the new axis-1 and axis-2 become 

ATf cos2 (j) + AT| sin^ - sin (/) cos (j){Ti, T2) 
(T3)2 + ((Ti)sin0+(T2)cos0)2 ' 
2_., ATf sin^ + ATI cos^ + sin cos 0(Ti, 

(T3)2 + ((Ti)cos0-(T2)sin0)2 ' 

where (T,,T,-) = {TiT, + T,T,) - 2(T,)(T,-). 



4. Results 



4.1. Numerical method 

The mean-field Bose-Hubbard Hamiltonian in ([3]) has been used to examine the phase 
transition between the superfiuid and Mott-insulator phases [S], with ip^. denoting the 
order parameter for the a mode. The superfiuid phase for the a component is identified 
with iptj 7^ 0. In the superfiuid state the tunnelling term Jo- is large and dominates 
the Hamiltonian. As a result the ground state corresponds to the single particle wave 
function of all cr-type atoms extended over the whole lattice, with each site being a 
coherent superposition of Fock number states [3]. In the Mott phase, on the other 
hand, the interaction term dominates so that the ground state exhibits minimal number 
fiuctuation and corresponds to a product of atom Fock number states at each lattice 
site, which in turn gives ipa = ^ l3j. 

We have performed numerical diagonalization of the mean-field Hamiltonian ([3]) by 
using a set of states expanded in terms of the product of individual atom number states 

N N 

\^) = C„onAl^o)|^A)- (9) 

no=0 nA=0 

While performing this diagonalization, two different regimes with respect to the same 
parameter P must be carried out. One is for a positive antisymmetric coupling, with 
a corresponding antiferromagnetic ground state, where individual spins are anti-aligned 
due to spin-exchange interaction. The other case is ferromagnetic for a negative spin 
exchange interaction. In addition, we explore the dependence of our results on the 
small, but non-vanishing lattice parameter 6, which introduces a spin dependent lattice 
potential. 

We study the parameter regions corresponding to those considered in Ref. [8]. 
The values of the parameters in Hamiltonian (I3l), which are needed for numerical 
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computation, are thus read from the Fig. 1 of Ref. [5], with J/U picked to ensure 
the system have full access to the n = 2 Mott regime, but barely enters the n = 3 Mott 
phase. 6 is taken to be small and 6 values used are for the range of < ^ < 1. We 
study the degenerate (6* = 0) and non-degenerate cases (6* 7^ 0) separately. From the 
initial values of the order parameters ipo and ipA we compute the diagonal basis and 
the corresponding ground state. This ground state then allows us to calculate the new 
order parameters and to compare with the initial values. This procedure is iterated to 
reach a self-consistent solution, with which it becomes straightforward to calculate the 
expectation values and the second moments of the operators in 

To conveniently calculate the squeezing parameter (Cj), we use the average total 
occupation number (n) for each type of interactions to label the different phases instead 
of relying on the total number of atoms N (per site). This implicitly assumes that 
the squeezing parameter (I7j) remains a valid criterion of quantum entanglement even 
for non-integer occupation numbers such as in the superfiuid phase. This assumption 
does not introduce any inconvenience in a Mott phase since the ground state consists 
of Fock states with equal total number of particles, i.e., definite spin and thus (h) 
becomes an integer. In the superfiuid phase, we justify the use of a non- integer (n) 
in the following manner. In this section, we calculate the squeezing parameter in two 
different ways for each case. The first method uses (n) directly for the entanglement 
measure. The second method is analogous in form, but only uses integer values of (n). 
For the superfiuid phase, instead of talking about separability for states with different 
total number of particles, we focus on the subspace no + n\ = n block and investigate its 
correlation. This becomes a meaningful measure when the block we use is the one with 
the nearest integer total number of particles to (n). This method has a similar nature 
as the superselection rules mentioned in Ref. [18] and in Ref. [30] since the projection of 
the Hilbert space onto a subspace of fixed particle number is considered. Both methods 
are found to give similar behaviors for the superfiuid and the Mott insulator phases. 
We provide results from the first method in our discussion because they respect the 
collective nature of the superfiuid state and emphasize particle number fluctuations. 

There also exist states for which spin squeezing parameter cannot be readily used 
to characterize their correlation properties. An example is the maximally entangled 
states (MES) in Ref. [31], which are not squeezed spin states according to the criterion 
in ([7]). In this case, it is inadequate to talk about squeezing, since the uncertainty 
in the perpendicular components to the mean isospin vector are meaningless as the 
denominator for the squeezing measure ([7]) vanishes for all axes. In addition, there exist 
other states, although whose averaged mean isospin are nonzero, the expectation values 
for the two components in the denominator might vanish, also making the spin squeezing 
parameter not well defined. In our studies, we find that these states happen only 
in certain Mott phases, where exact wave functions are available either analytically in 
the spin [23] or Fock basis [8]. As such, their quantum entanglement properties can be 
discussed directly using other criteria. 

In order to quantify the pairwise quantum correlations both in the superfiuid and 
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Mott-insulator regimes, in addition to the squeezing parameter, we use the well-known 
criterion called concurrence [28]. For a given two-party state p, this measure is equal to 

=max{0,Ai-A2-A3-A4}, (10) 

where Aj's are the square roots of the eigenvalues of pp in decreasing order where 

p= {ay^ay)p*{ay®ay). (11) 

For the n = 1 Mott-insulator phase, this measure is trivial since there is only 
one particle present. When it comes to the n = 2 Mott phase, concurrence clearly 
quantifies pairwise correlations between the two atoms at the same lattice site. In the 
superfluid phase, the ground state is a superposition of Fock states with different number 
of atoms or isospin states with different isospins, we again focus on the subspace with the 
nearest integer total number of particles uq + = n. If the nearest integer is smaller 
than two, then concurrence is zero. If it is equal to two, the concurrence is simply 
calculated. When it is equal to three, the three-particle ground state is symmetrized in 
the first quantization picture and we use reduced two-body density matrix to calculate 
concurrence. 

We report below our investigation of quantum entanglement in our model system 
for the two regimes: antiferromagnetic and ferromagnetic interactions. 



4.2. Ferromagnetic regime 

For ferromagnetic interaction with P < 0, for ^ = 0, and a fixed J/U value, the 
dependence of the order parameters ipQ and ip;^ on the quantity p/U is shown in Fig. [TJ 
We determine the phase of the system for any p/U value by looking at the order 
parameter of each component. 

li . . . . . 1 




Figure 1. The dependence of the order parameters on the value of fi/U for 9 = 0, 
J/U = 0.455 X 10~^, and P/U = —0.926 x 10~^ in the ferromagnetic regime. The 
vanishing of the order parameters matches closely with the appearance of Mott- 
insulator phases for the corresponding component. 
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Quantum correlations on a single lattice site is evidenced by evaluating the 
squeezing parameter ([7]). In the superfluid regime, numerical calculations, taking into 
account the minimization with respect to coordinate rotations, yield ^.f > 1- So there 
is no particle entanglement in the superfluid phase of a ferromagnetically interacting 
system when ^ = 0. Although, this situation deserves to be more carefully analyzed for 
those values of fi/U that correspond to Mott-insulator phases. 

The spin squeezing parameter is not defined for the zero particle (n = 0) ground 
state |0,0), the trivial case of no particle entanglement without any particles. When 
6 = 0, the single particle (n = 1) ground states |10) and |01) are degenerate [8] and can 
be written as \g) = cosx|01) + sina;exp (i?/)|10), where x,y E [0, 27r] are arbitrary angles, 
parameterizing the manifold of the ground state family. We find (Ti) = (1/2) sin 2x cos y, 
(T2) = (—1/2) sin 2x sin?/, and (T3) = (—1/2) cos2x. The spin fluctuations are (ATf) = 
(1/4)(1 -sin2 2xcos2?/), (Ar|) = (1/4)(1 - sin^ 2x sin^ y), and (AT|) = (l/4)(sin2 2x). 
Thus we obtain C,f = 1 in any direction i = 1,2, 3, for any member of the ground 
state manifold. The ground state, expressed in the spin representation [21], could be 
written as an arbitrary superposition of |T = 1/2, Ti = ±1/2) spin states. We write 
\g) = \x,y) = cos (x/2)|l/2, 1/2) + sin (x/2) exp (i?/)|l/2, — 1/2) for the ground state 
in spin representation. Projection of the total spin onto the (x, y) direction gives the 
spin component S^^y = sinxcos?/T2 + sina:sin?/T3 + cosxTi, whose eigenstate is \x,y) 
with eigenvalue 1/2, such that Sx,y\x,y) = {l/2)\x,y). Such a state is called a coherent 
spin state (CSS) [T9]. The ground state \g) is identified as a pure state of a spin-1/2 
system, and as such is a CSS. There exists no other spin to be correlated with, so that 
\g) cannot be a squeezed spin state (SSS). Particles in a CSS are correlated as all spin 
1/2 constituents atoms are pointing along the same direction; although they remain 
separable, i.e., they are not entangled. 

On the other hand, the n = 1 Mott state could become mode entangled [32] for 
some a and /3. Mode entanglement is a different concept from particle entanglement 
considered here and could be useful for different applications [32]. It corresponds 
to entanglement in the second quantization picture, while particle entanglement is 
associated with the inseparability of the wave function, or density matrix, in the first 
quantization. 

Similarly, the ground states for the n = 2 Mott phase are also degenerate for 6 = 0. 
As such they form a manifold represented by \g) = cosx|ll) + sinxexp {iy)\b), where 
\b) = (|02) + |20))/v^. In this case, (Ti) = sin2xcos?/ and (T2,3) = 0. The variances 
are calculated to be (ATf) = 1 — sin^ 2x cos^ y, {AT2) = cos^x, and (AT|) = sin^a;. 
becomes either undetermined (a 0/0 form) or 00 due to vanishing denominators. If we 
calculate after a coordinate rotation by about the axis-3, we find ^f,. Minimizing it 
with respect to 0, we finally get (^i/)min = 1/(2 sin^ x cos^ j/) with its minimum value at 
(j) = ±71 / 2. We find ^| = 1/(2 cos^ x cos^ y) and ^| = 1/(2 sin^ x cos^ y). For some values 
of X and y, < 1 is satisfied. Hence, particle entanglement exists for some members of 
the ground state manifold. This is consistent with the fact that each degenerate ground 
state 1 11) and \b) is particle entangled. For parameters x and y specifying a dominant 
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contribution from a particular degenerate component in \g), particle entanglement is 
expected. In the spin representation, the ground state is an arbitrary superposition of 
\T = l,Ti = ±1,0). In contrast to the spin-1/2 case of the n = 1 Mott phase, now 
SSS (squeezed spin state), where all particles are entangled, can be found in the ground 
state family. 

When we analyze ferromagnetic regime by calculating the concurrence in light of 
the discussion in Sec. 14. ![ it is found to be zero for all fi/U values except those for the 
n = 2 Mott phase. In this case, the ground state is an arbitrary superposition of two 
degenerate maximally entangled states, with the concurrence for each state being equal 
to one. But the concurrence for the ground state manifold mentioned above becomes 
C{\g)) = [1 — (1/2) sin^(2x) cos(2y)]^/^, which is larger than zero for some values of x 
and y. And this indicates the possibility of pairwise entanglement for certain ground 
states. 

Now we look at the situation when 6 takes a small but nonzero value. In this case 
the relations Jq ~ = J, Uq ^ U\ = U, and K = U + P remain valid. However, 
the parameter 6 becomes nonzero, due to the splitting between the two ground state 
modes: a = and a = A. This causes the dependence of the order parameters on fi/U 
to change as illustrated in Fig. [2l Note the difference between the order parameters for 
the two modes. 

i.2r • • • • • 1 




0.5 1 1,5. 2 2.5 3 



Figure 2. The dependence of order parameters for the two modes vs ///[/ for a small 
nonzero 9 in the ferromagnetic regime with J/U = 0.625 x 10~^, P/U ~ —0.926 x 10~^, 
and 6/U — 0.327 x 10^^. The solid line denotes ipA while the dashed line refers to ipQ. 

For the general case with 6^0, performing minimization over the axis rotations 
shows that the optimal squeezing occurs for the unrotated coordinate axes. By 
examining the spin squeezing parameter as a function of fi/U numerically, we find 
that particles are not entangled in the superfluid regime. We thus look for particle 
entanglement in the Mott phases. 

With a small 6, the degeneracy in the ground states in the Mott phase is removed. 
In the n = 1 Mott-phase, the ground state becomes \g) = \no = l,nx = 0). For this 
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state, the mean spin is along the direction of the axis-3 with (T3) = 1/2 and T12 = 0. 
The spin fluctuations are given by (AT^g) = 1/4 and (ATg^) = 0. Employing a rotation 
by about the axis-3, we find that ^f, 3/ = 1, i.e., the ground state is a CSS. 

For the n = 2 Mott-phase, we have a non-degenerate ground state of the form 
\g) = a|02) + 6|20) [8] with 



(12) 




(13) 



For such a state, as in the n = 1 Mott phase, the mean spin is pointed along the axis-3 
with (Ti^2) = and (T3) = 6^ _ Their corresponding fluctuations are found to be 

To determine 



(AT2) 



6)2/2, (AT|) = (a - 6)2/2, and (AT32) = 1 - (^2 



,2\2 



the optimum noise reduction and spin squeezing, we minimize over rotations about the 
mean spin (axis-3) direction by an angle 0. It is sufficient to consider either one of the 
rotated 1' or 2' axes so that a single rotation angle dependent spin squeezing parameter 
^2 can be found as 

1 + 2abcos2(j) 



(14) 



(62 - a2)2 ■ 

Its minimum occurs at = ±7r/2 such that ^^^/2 ~ ~ 2ah)/{}P' — a^Y. Assuming a 
small 5/P, we find ^^^/2 ~ + 0(((5/P)2), in agreement with numerical calculation 
reported in Fig. O Thus, the ground state is particle entangled and spin squeezed. 

We again calculate the concurrence values for the phases under consideration. It 
becomes zero everywhere except n = 2 Mott phase. In this situation C{\g)) = 2\ab\ 
and for small 6/P values C{\g)) ~ 1 — 0((5/P)2). So that the results are in complete 
agreement with those of squeezing parameter. 



4.3. Antiferromagnetic regime 

In this case, the atomic interaction parameter P is positive. In Fig. [3l the order 
parameters are plotted as a function of fi/U at 9 = 0. 

Similar to the ferromagnetic regime, we first test the existence of spin squeezing for 
6 = 0. The corresponding minimum squeezing parameter, ^| for the fixed axes is shown 
in Fig. m In contrast to the ferromagnetic case, squeezing is observed for the superfiuid 
phase as well. In numerical calculations, we also rotate the coordinate system to see 
whether correlations can be enhanced for some angles. The optimum squeezing is found 
to occur for the fixed axes configurations. 

In the n = 1 Mott-phase, the ground state is a coherent superposition of |10) and 
|01), which identifies a manifold of any pure state for spin-1/2. The only difference 
being the quantization axis, it lies along the axis-2, instead of the axis-1. Hence our 
conclusions for the ferromagnetic case remain applicable. The ground state family is a 
general CSS and exhibits no squeezing, although mode entanglement can be present. 
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Figure 3. The dependence of the order parameters on the values of fi/U for 9 = 
in the antiferromagnetic regime with J/U = 0.455 x 10~^ and P/U = 0.926 x 10~^. 
The nonzero valued order parameters indicate superfluid phases for the corresponding 
components. 




CM CM 
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0.5 1 1.5 2 2.5 3 

U/U 

Figure 4. The minimum squeezing parameter ^| for the fixed axes configuration in the 
antiferromagnetic regime with 6 = 0, J/U = 0.455 x 10~^, and P/U = 0.926 x 10~^. 
^2 < 1 denotes spin squeezing for the axis-2. 

The n — 2 Mott insulator state in the antiferromagnetic case, however, is 
significantly different from the ferromagnetic case considered earlier. It is no longer 
degenerate as before, and becomes uniquely determined as 

|5) = ^(|20) + |02)), (15) 

instead. For this special superposition state, the mean isospin vector becomes zero, with 
(^1,2,3) = 0. Spin fluctuations are found to be (ATpg) = 1 and {AT2) = 0. Given in 
the second quantization form and in the occupation number representation, the mean 
number of particles in each mode (0, A) is 1 and the state is mode entangled. In the 
first quantization, denoting single particle wave functions as ^'i^ for particles i = 1,2 
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in modes a = 0, A, \g) is found to become \g) = (l/-\/2)(\l/io^I/2o + ^ia^2a)- This state 
has maximum quantum correlation among the particles and can be identified as a MES 

m- 

In order to compare the results measured in terms of the calculated concurrence, 
we show in Fig. [5] the dependence of concurrence as a function of fi/U. 



1 
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Figure 5. Concurrence vs fi/U in tiie antiferromagnetic regime with 9 = 0. The 
presence of pairwise entanglement is assured if the value of concurrence becomes larger 
than zero. 



The presence of particle entanglement in the superfluid phase is reflected by the 
nonzero values of concurrence for the corresponding fi/U values as shown in Fig. [5l 
Having a concurrence of one in the n = 2 Mott phase corresponds to the presence of a 
maximally entangled ground state. 

As is done previously for the ferromagnetic small nonzero 9 value can be 

introduced and the system parameters are changed accordingly. The corresponding 
graph for the order parameters as functions of /i/f/ are shown in Fig. [6l 

Following the earlier procedure, the minimum squeezing parameter ^| is also plotted 
against /i/f/, with the optimized values, corresponding to the flxed coordinate system 
shown in Fig. [71 

As in the case of = 0, spin squeezing is found to exist for the superfluid phase 
almost with the same strength. On the other hand, although spin squeezing is detected 
in the n = 2 Mott phase, it is reduced with a nonzero 6. The corresponding ground 
state for the n = 2 Mott phase is the same as in the ferromagnetic case. The MES of 
the ^ = case for the antiferromagnetic interaction becomes a partially entangled state 
when a small nonzero 6 is introduced. 

The results from the calculated concurrence as shown in Fig. [8] are in complete 
agreement with those from the squeezing parameter. Squeezing is present in the 
superfluid phase and the maximal entanglement in the n = 2 Mott phase becomes 
partially entangled with the introduction of a small but nonzero 6. 
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Figure 6. The order parameter tpo for the antiferromagnetic regime at a small 6 with 
J/U = 0.455 X 10-\ P/U = 0.926 x lO'^, and S/U ^ 0.327 x lO'^. V'a = for these 
values of interaction parameters. 



1.1 




Figure 7. The minimum squeezing parameter ^| in the antiferromagnetic regime at a 
small in the fixed coordinate system with J/U = 0.455 x fO"\ P/U 0.926 x fO"^^ 
and 6/U = 0.327 x 10^^. 



5. Conclusion 

In summary, we have investigated quantum correlations between spin-1 bosons with 
coupled ground states in optical lattices. Both ferromagnetic and antiferromagnetic 
interactions are considered based on a model, initially developed in Ref. [8], that 
we believe can be readily adopted to current experimental systems. In addition to 
characterizing quantum correlations in various quantum phases in terms of coherent 
and squeezed spin states, and addressing both particle and mode entanglement, the role 
of lattice parameter in the familiar lin-^-lin configuration is examined. 

We have shown that for ferromagnetic interactions isospin squeezing (or multi- 
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Figure 8. Concurrence at a small but nonzero 9 in the antiferromagnetic regime. It 
is seen that part of the superfluid phase contains entangled particles. 



particle entanglement) is absent in the lattice model of spin-1 bosons in the superfluid 
phase. The one particle Mott phase is in fact in a CSS, which is not particle entangled, 
although it displays significant mode entanglement. The two particle Mott state may 
contain SSS and entangled particles, if one of the degenerate component in the ground 
state manifold is made to dominant. It can be steered into a particle entangled state 
by introducing a nonzero 6 to lift the degeneracy, while the CSS of the n = 1 Mott 
phase or the superfluid phase remains unentangled. The path to quantum entanglement 
is through the well known single axis twisting type nonlinear interaction [T9] for the 
degenerate {6 = 0) case. With a nonzero 6, quantum entanglement is generated from 
a generalized LMG interaction, which includes a two-axis twisting type of spin-spin 
nonlinear interaction. 

For antiferromagnetic interactions, spin squeezing and particle entanglement is 
found in both the n = 2 Mott and superfluid phases. In the n = 2 Mott state we 
find maximally entangled particles. Introducing a nonzero 6 reduces this to a partially 
entangled state, and thus decreases particle correlations. 

We compared the results of the squeezing parameter ([7j) with those of the 
concurrence (11 01) for each type of interaction and lattice configuration. They are in 
complete agreement in demonstrating the presence or absence of entanglement for the 
different phases. 

For the system under consideration, we have investigated the potential 
ground states and the corresponding quantum correlations via examining entangle- 
ment/squeezing properties. Depending on the interaction parameters of the system, 
abrupt changes may occur if one considers the behavior of entanglement properties. 
One can introduce symmetry breaking perturbations to the Hamiltonian ([5]) to remove 
the degeneracy present in the various ground states. This can be done via including mag- 
netic fields and Raman pulses with which adjustments to the ground state populations 
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in any particular spin components can be made [33] • As a specific example, generation of 
a coherent superposition of degenerate states (in this case Zeeman sublevels Mp = ±1) 
by stimulated Raman adiabatic passage scheme is demonstrated experimentally in Ref. 

m- 
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